Particle production in central Pb+Pb collisions at 158 A GeV/c 
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The NA44 experiment has measured single particle inclusive spectra for charged pions, kaons, and 
protons as a function of transverse mass near mid-rapidity in 158 A GeV/c Pb+Pb collisions. From 
the particle mass dependence of the observed tut distributions, we are able to deduce a value of 
about 120 MeV for the temperature at thermal freeze-out. From the observed ratios of the rapidity 
densities, we hnd values of the chemical potentials for light and strange quarks and a chemical 
freeze-out temperature of approximately 140 MeV. 

PACS numbers: 25.75.-q,13.85.-t,25.40.Ve 



I. INTRODUCTION 



The goal of heavy-ion collision research at ultra- 
relativistic energies is to understand the fundamental 
properties of nuclear matter under extreme conditions of 
high particle and energy density. When these densities 
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and temperatures in heavy-ion collisions are sufficiently 
high, quarks should no longer be confined in hadrons 
but instead move freely through the entire volume. The 
colliding system expands and cools, driven by the large 
internal pressure generated by the strong interaction of 
the constituents. The quarks eventually become confined 
again, forming hadrons that can be detected by experi- 
ments. The hadron yields and momentum distributions 
provide information about the temperature, density, and 
dynamics of the later stages of the collision. From this 
inferences about the conditions at earlier stages can be 
made. 

The multiplicity of the produced hadrons ranges from 
several hundred to a few thousand in relativistic heavy- 
ion collisions at ^/snn ~20 GeV, and their rapidity den- 
sity cannot be explained by simple superposition of p+p 
collisions, as final state rescattering may play an impor- 
tant role. The transverse momentum distributions of the 
hadrons display a Boltzmann-like shape, inspiring a ther- 
modynamical interpretation of the momentum distribu- 
tions. 

In this report, we present invariant cross sections of 
tt^, K , p, and p, as a function of transverse mass near 
mid-rapidity in central 158 A GeV/c Pb+Pb collisions. 
NA44 has previously reported the particle yields for pi- 
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FIG. 1: The schematic view of NA44 spectrometer. 



ons, kaons, and protons in p+A to S+A collisions [|, 
kaon production in Pb+Pb Q, and antiproton and an- 
tideuteron production in Pb+Pb J5|. Because of the 
excellent particle identification and systematic measure- 
ments from small to large systems, NA44 is well-suited 
to study the thermal and chemical properties of the 
hadronic system at freeze-out. The dependence of the 
thermal parameters on the colliding system size will also 
be discussed. 



II. NA44 EXPERIMENT 

The NA44 focusing spectrometer is designed to mea- 
sure one and two particle momentum distributions of 
charged hadrons near mid-rapidity Q . A schematic view 
of the spectrometer is shown in Fig. [l| The target is a 
lead disk of 2 mm thickness (3.4% interaction probabil- 
ity for a lead beam). A Cherenkov gas beam counter 
(CX) and a beam veto counter (CX-veto) are placed 
upstream of the target. CX selects single beam parti- 
cles (rejecting events with multiple simultaneous beam 
particles) up to a rate of 2 x 10 6 ions/second while the 
CX-veto signals the presence of any beam halo. A scin- 
tillation multiplicity counter (TO) is 10 mm downstream 
of the target, and produces a pulse height proportional 
to the number of charged particles impinging upon it. 
TO covers 20% of the azimuthal angle and has a pseudo- 
rapidity range of 0.6<?7<3.3. This contains mid-rapidity 



(y-mid =2-9) for the Pb+Pb collision. TO is used to select 
collision centrality at the trigger level. 

Two dipole magnets (Dl and D2) analyze the particle 
momenta. Three super conducting quadrupole magnets 
(Ql, Q2, and Q3) focus the particles. The PID/tracking 
detectors measure hit positions after the magnetic anal- 
ysis, and consist of strip chambers and scintillation ho- 
doscopes, along with two threshold gas Cherenkov coun- 
ters, CI and C2. Three hodoscopes, H2, H3, and H4 
consist of 60, 50, and 60 scintillation counters, respec- 
tively. Only H2 and H3 are used in this analysis. The ho- 
doscopes provide position information for tracking, with 
vertical hit information determined from photomultiplier 
tubes on each end of every counter. The hodoscopes 
also measure time-of- flight for hadron identification [0. 
One pad chamber (PC) and two strip chambers (SCI and 
SC2) provide 0.3 mm position resolution for track recon- 
struction. CI and C2 contain Freon-12 at 2.7 atm and 
N2 at 1.3 atm respectively, and identify particles at the 
trigger level. The trigger requires a valid beam, a high 
multiplicity detected in TO (central collision) , at least one 
track in each hodoscope, and the absence of an electron 
(all runs) or pion (for K and p runs). 

The timing signal from CX is used for the time-of-fiight 
(TOF) start with a resolution of 35 ps. The average of 
the top and bottom PMT of the hodoscope is used for the 
stop signal. The overall TOF resolution is approximately 
100 ps and 80 ps for H2 and H3, respectively. 

Two magnetic field settings of the spectrometer 
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FIG. 2: The 7r, K, and p acceptance in rapidity and pr- In this analysis, pions at the 8 GeV/c 44 mrad angle setting are 
selected with pt>0.35 GeV/c. 
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FIG. 3: An example of particle identification in the 4 GeV/c setting. The contour plot shows pulse heights of CI and C2 
(left). The region (a) is electron veto and pion required by CI and C2 pulse heights. The region (b) is electron and pion veto 
by CI and C2. The right-hand plots are plots of the time-of-flight vs. track momentum for events inside the boxes labeled (a) 
and (b) in the left-hand plot. The pair of plots on the right are TOF at H3 as a function of momentum with CI and C2 cuts. 
The curves show ±2.5 sigma of TOF resolution from the peak. Particles between the curves are selected. 



are used, yielding nominal momentum of 4 GeV/c 
(3.2<p<5.2 GeV/c) and 8 GeV/c (6.0<p<9.8 GeV/c). 
Two angular settings, 44 mrad and 129 mrad are used to 
give a large range in transverse momentum {pr)- Com- 
bining these settings yields transverse momenta from 
p T = to 0.9 GeV/c (0 to 1.6 GeV/c) for 4 GeV/c (8 
GeV/c) setting. Figure || shows the NA44 acceptance in 
Pt and rapidity for pions, kaons, and protons. The NA44 
acceptance includes mid-rapidity (y mid =2.9) for Pb+Pb 
collisions at 158 GeV/c per nucleon. 

The invariant cross section is determined by the fol- 
lowing analysis procedures: track and momentum recon- 
struction, particle identification of tracks, a selection of 
high multiplicity events, an acceptance correction via a 
Monte-Carlo simulation, normalization by summing of 
beam scalers, and multiplicative correction factors for 
cuts. 

The hit positions on PC, H2, H3, SCI, and SC2 were 
used to reconstruct tracks by fitting them with a straight 
line using a ^-minimization method. As is clear from 



TABLE I: The number of particles selected by the cuts for 
each spectrometer setting. The event fraction is top 3.7 ± 
0.1%. 
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Fig. [j], the tracking detectors lie outside the magnetic 
field. The momentum resolution Sp/p is 0.2% for all par- 
ticle species. 

Particle identification uses a combination of pulse 
heights in CI and C2 with time-of-flight from CX to 
the hodoscopes. The Cherenkov threshold momenta of 
C1(C2) for e, tt, K, and p are 0.014, 2.5, 8.9, and 16.9 
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FIG. 4: The pulse height distribution of the TO multiplicity 
counter for valid beam events (1) and events that satisfy the 
TO centrality trigger (2). The shaded region is used for the 
analysis. 



(0.020, 5.3, 18.3, and 34.8) GeV/c, respectively. In the 4 
GeV/c spectrometer setting, pions are separated from K 
and p by CI and C2, and then K/p separation is done 
by TOF. In the 8 GeV/c setting, e and it are vetoed by 
C2, and then K/p separation is done by CI. The e/it 
separation is not able to be done by CI and C2 in the 
8GeV/c setting, however the contamination from elec- 
trons is negligible for pr>0.35 GeV/c. The statistics for 
each particle are shown in Table ffl. 

Because it has a longer flight path in addition to its 
better timing resolution, H3 rather than H2 provides the 
primary TOF information. Figure [| shows an example of 
particle identification in the 4 GeV/c setting. The left- 
hand side shows a scatter plot of pulse heights in CI and 
C2. Events containing an electron fire C2 - these events 
are rejected. The events selected by the box labeled (a) 
and (b) on the left-hand side of figure are shown on plots 
of TOF as a function of momentum on the right-hand side 
of the figure. Particles between the curves are selected. 
The contamination for each particle species is less than 
3% in both momentum settings. 

We trigger on central collisions using the interaction 
counter, TO, which measures the total amount of light 
produced by charged particles crossing its scintillators. 
Figure |^ shows the pulse height distribution of the TO 
for valid beam events (1) and central triggers (2). The 
shaded region is used for the analysis. The width of the 
valid beam distribution comes from the production of 
delta electrons. Since TO covers only 20% of the az- 
imuthal angle, it is not possible to count all low mul- 
tiplicity Pb+Pb interactions and we cannot measure the 
Pb+Pb interaction cross section. Therefore we scale the 
measured p+Pb cross section by a geometry factor: 
this fraction is denoted by A. The value of A is ob- 
tained by scaling the interaction probability for p+Pb 
collisions (the nuclear interaction length) to the interac- 



tion probability for Pb+Pb collisions by a geometry fac- 
tor of (2083 +2083) 2 /(l* +2083) 2 , i.e., A=2.93. The nu- 
clear interaction length of Pb is 194 g/cm 2 and density of 
Pb is 11.35 g/cm 3 . Therefore, the interaction probability 
of p+Vh of 2 mm Pb target is 0.2/(194/11.35) = 0.0117. 
The total Pb interaction probability is 3.4%(= AxO.0117) 
for 2 mm Pb target. 

We then find that the shaded region in Fig. || corre- 
sponds to 3.7+0.1% (statistical only) of the Pb+Pb inter- 
action cross section. This region corresponds to a region 
of impact parameter (6)~<3 (RMS~1) fm, estimated by 
applying the TO acceptance cut to the multiplicity dis- 
tribution from RQMD (v2.3) 0. We note that the TO 
covers only 20% of the azimuthal angle, therefore the 
event centrality selected by TO has a weak sensitivity to 
multiplicity, that is, to impact parameter. 



III. ANALYSIS 

In order to construct the invariant cross section, the 
effect of the acceptance is estimated by a Monte-Carlo 
(MC) simulation based on TURTLE Jn|, which passes 
generated tracks through the detectors. The acceptance, 
reconstruction efficiencies, momentum resolution, and 
loss of pions and kaons due to in-flight-decays are esti- 
mated by comparing generated tracks with reconstructed 
tracks. The momentum spectrum and rapidity distribu- 
tion for each particle species in Pb+Pb collisions is taken 
from RQMD (v2.3). The Monte-Carlo events were ana- 
lyzed with the tracking programs exactly as applied to 
the physics events. The corrections are given by the ra- 
tio of the generated yield to the reconstructed yield as a 
function of m-r inside the rapidity acceptance. The abso- 
lute cross section as a function of Tot in a rapidity range 
Ay is given by simplify this equation. 
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where n is the number of measured tracks, N event is the 
number of collisions, h mc is the number of particles as 
input of the MC simulation and h r is the number of tracks 
reconstructed in the simulation. 

Secondary hadronic interactions of protons and an- 
tiprotons in the spectrometer material does not distort 
the shape of the momentum distributions but reduces the 
yield [||. The correction factors are estimated to be 1.11 
for protons and 1.17 for antiprotons, respectively. 

The corrected momentum distribution is then normal- 
ized using the number of beam particles, the interac- 
tion probability, the fraction of interactions satisfying the 
trigger, the measured live time of the data-acquisition 
system, and losses due to particle identification cuts. 

Separation of protons and kaons from pions is done via 
Cherenkov veto trigger, but in the process some protons 
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and kaons are rejected in events that also include a pion. 
The veto factor (i.e., the correction for this loss) for each 
particle and momentum setting is estimated by compar- 
ing the number of the particles with no Chcrcnkov veto to 
the number of particles with Cherenkov veto, using real 
data without the Cherenkov veto trigger. Since the veto 
factor depends on relative particle ratios in spectrometer, 
it must be independently estimated for each momentum 
and angle setting. 

The transverse mass distribution is obtained from the 
two angle settings (44 mrad and 129 mrad). The data 
at 129 mrad have smaller (typically one half to one 
third) systematic errors in the absolute normalization 
than those at 44 mrad. Therefore, in this analysis, the 
data at 44 mrad were scaled to match those at 129 mrad 
in the mj region where they overlap. The distributions 
from both settings were fit by exponentials with the same 
inverse slope parameter. The dN/dy is the sum of the 
normalized cross section in each measured mj- bin within 
the acceptance plus an extrapolation that uses the fitted 
exponential distribution beyond the measurement region. 
This extrapolation contributes less than 1% of the yield 
for pions but up to 7% for protons with the 4 GeV/c 
setting. 



IV. SYSTEMATIC ERRORS 

The inverse slope parameters are extracted by fitting 
the experimental distribution with a single mj exponen- 
tial function (see Eq. ([!]) in the section VI). Systematic 
errors on the inverse slope parameters are estimated by 
dividing the data into two px ranges for the fits. The 
systematic errors are shown with the inverse slope pa- 
rameters in Table Q. 

The absolute normalization procedure results in sig- 
nificant uncertainties in dN/dy. Contributions include 
uncertainties on the following; the data-acquisition live 
time (5%), the stability of the TO pulse height used to 
select the event fraction (2%), the Cherenkov veto fac- 
tor (<5%), the correction for particle loss due to quality 
cuts (<5%), and the extrapolation beyond the measure- 
ment region (<1~7%). The last three vary with particle 
species and spectrometer setting. We assume that the er- 
rors are uncorrelated and so add them in quadrature to 
get the total systematic error. The errors for each parti- 
cle species and setting are shown with results in Table ||. 
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FIG. 5: The correction factors for proton and antiproton as a 
function of mr— mass. The solid (open) circle corresponds to 
44 (129) mrad setting. The horizontal bars above and below 
each point show the systematic errors from uncertainties in 
the hyperon/proton ratios. 



E + and £°. Since the E° lifetime is short, its yield was 
included in the A yield. The dN/dy for various parti- 
cles near mid-rapidity in central collisions was taken from 
RQMD(v2.3). 

The correction factor as a function of mj is 



C p {m T ) 
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where N Paii is the number of the original protons and 
N PA and N Ps represent the protons from A decays, and 
S + decays that track through the spectrometer. Figure || 
shows the correction factors as a function of m^— mass. 
Since these depend on the ratio of A and E + to pro- 
ton, the production ratios were varied. The circles are 
the correction factor estimated from A/p and S/p ratios 
of RQMD, while the lines assume yield ratios larger or 
smaller by a factor 1.5. 

The correction procedure contributes to the systematic 
uncertainties on the inverse slopes of proton and antipro- 
ton by about 7% for the 4 GeV/c setting and about 3% 
for the 8 GeV/c setting. The uncertainties due to feed- 
down correction on dN/dy are approximately 15% for 
protons and 19% for antiprotons in both 4 GeV/c and 8 
GeV/c settings. 



V. FEED-DOWN FROM WEAK DECAYS 

NA44 has no vertex detectors in the target region, so 
decay protons from some fraction of the short-lived par- 
ticles such as strange baryons will be detected as protons 
in the spectrometer. This contamination of protons and 
antiprotons is not negligible. 

The same Monte-Carlo used for acceptance corrections 
was also used to estimate the weak decay yields from A, 



VI. RESULTS 

The invariant cross sections measured by the NA44 
spectrometer in Pb+Pb collisions are shown in Fig. ||. 
They are plotted as a function of tot in the specified 
rapidity regions. The tot distribution is generally well 
described at these energies by 
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FIG. 6: Invariant cross section of n, K, and p in Pb+Pb collisions as a function of mr — mass from the 4 GeV/c and 8 GeV/c 
setting. The open (solid) marker shows positive (negative) charged particle. The rapidity range for each particle is shown in 
the figure. The proton and antiproton distributions are corrected for feed-down. The horizontal bars above and below each 
point show systematic errors from the feed-down correction and correspond to the horizontal bars in Fig. Note that the 
proton and antiproton are scaled by a factor of 0.1 in the vertical direction. 



TABLE II: The inverse slope parameters and dN/dy with statistical errors (first) and systematic errors (second). The pion 
dN/dy is obtained from the data of 4 GeV/c and 8 GeV/c setting at both 44 mrad and 129 mrad. The values for proton and 
antiproton including the feed-down correction are shown as p c and p c , and systematic errors include errors from feed-down 
correction. 
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7T + 
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207 ± 3±| 
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126 ± 2±f 


42.3/13 






2.4-4.1 








150 ± 3 ± 8 


7T~ 


2.4-3.1 


0.28-1.20 


201 ± 3±? 


63.0/25 






3.1-4.1 


0.00-0.56 


109 ± 2±|° 


67.9/13 






2.4-4.1 








162 ± 2 ± 10 


K+ 


2.3-2.9 


0.00-0.32 


221 ± 9±2° 


17.1/9 


24.0 ±0.2 ±2.1 




2.4-3.5 


0.00-0.84 


246 ± 2±f 


79.1/28 


27.5 ±0.2 ±1.6 


K~ 


2.3-2.9 


0.00-0.32 


224 ± 7 ± 1 


7.7/9 


14.8 ±0.5 ±1.1 




2.4-3.5 


0.00-0.84 


228 ± 2±? 


49.2/28 


15.4 ±0.5 ±1.0 


V 


1.9-2.35 


0.00-0.28 


319 ±23±i 8 


10.4/6 


33.3 ±0.8 ±4.0 




2.3-2.9 


0.00-0.68 


296 ± 5±7 7 


74.1/16 


34.7 ± 1.0 ±2.2 


P 


1.9-2.35 


0.00-0.28 


303 ± 35±? 4 


9.0/6 


1.74 ±0.01 ±0.33 




2.3-2.9 


0.00-0.68 


300 ± 9±! 


18.5/16 


2.67 ±0.05 ±0.26 


Pc 


1.9-2.35 


0.00-0.28 


379 ± 33±| 8 


10.9/6 


27.3 ±0.5 ±5.2 




2.3-2.9 


0.00-0.68 


327 ±6± 2 ?, 


94.5/16 


28.5 ±0.2 ±4.6 


Pc 


1.9-2.35 


0.00-0.28 
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where A is a constant and T is the inverse slope param- 
eter. Table [fj lists the inverse slope parameters, the fit 
region, and the values of dN/dy. 

The slope parameters arc consistent with previous 
NA44 results (Ref. U)) within errors. The proton dN/dy 
is consistent with the values in Ref. . The statistics in 
our previous publication are 1/3 to 3/4 of the statistics 
used in this paper. We note that the NA44 collaboration 
published a paper on pion and kaon spectra focusing on 
"strangeness enhancement" Q . This analysis used a dif- 
ferent definition of the event fraction, based on the silicon 
multiplicity counter. The results are consistent within er- 
rors. 

The NA49 collaboration has reported dN/dy distribu- 
tions of 7r + , 7r~ , K + , K~ ,p, and p for the 5% most cen- 
tral Pb+Pb reactions Q [r| |l7| . These values are 
slightly higher than our data (top 3.7% event fraction). 
However, since they trigger on the energy of forward go- 
ing spectators their selection of low impact parameter is 
better than ours. The top 5% NA49 data corresponds to 
mean impact parameter (b}=2.2 fm [jlfi]]. On the other 
hand, the NA44 data corresponds to (6)~3 fm. Com- 
paring similar impact parameter regions, the values of 
dN/dy are consistent with the NA49 results within the 



We have obtained the invariant cross sections for pions, 
kaons, and protons in Pb+Pb collisions using data analy- 
sis procedures compatible with these used in Refs. ^ [|. 
Consequently, we can study the systematics of particle 
yields at CERN-SPS energy for many collision systems 
(450 GeV/c p+A, 200 A GeV/c S+A, and 158 A GeV/c 
Pb+Pb). 

Figure shows the inverse slope parameters for pions, 
kaons, and protons versus the dN c h/dy near central ra- 
pidity. All of the inverse slopes increase with dN c h/dy 
but this effect is more rapid for protons than for kaons 
and pions. The slopes of particles and antiparticles are 
equal except for protons have slightly larger inverse slopes 
than antiprotons for p+A and S+A collisions. While the 
uncertainties on the slope parameters for pions make a 
definitive statement difficult, the pion slope parameters 
in the region jtit— mass<0.6 GeV/c 2 in Pb+Pb are sim- 
ilar to, or less than, the slope in p+A and S+A collision 
systems. Perhaps, the lower values in the low px region 
for Pb+Pb reflect a larger enhancement at low pr- 

The kaon, proton and antiproton inverse slope param- 
eters increase with collision system size slowly, but more 
rapidly than the pion inverse slope parameters. The pro- 
ton inverse slopes are higher at mid-rapidity (y=2.3-2.9) 
than at more backward rapidities (y=1.9-2.3) in p+A 
and S+Pb collisions. 

Figure || shows the collision system dependence of 
charged pion, charged kaon, proton, and antiproton 
yields, and ratios of negative to positive charged particle. 
The pion rapidity density in Pb+Pb collisions increases 
by a factor of 71 ± 10 compared to p+Be collisions, and 
by a factor 5.1 ± 0.6 compared to S+S collisions. dN/dy 
for 7r~ and 7r + is compatible, within the errors for all 
collision systems. 

The mid-rapidity dN/dy of K+ (K~) in Pb+Pb colli- 
sions increases by a factor of 133+20 (105+10) compared 
to p+Be collisions and 6.4 ± 0.9 (6.2 ± 0.8) compared to 
S+S collisions. The system size dependence is stronger 
than for pions, i.e., the K/ir ratio increases from p+Be 
to Pb+Pb. 

Proton dN/dy increases with charged multiplicity at 
mid-rapidity more rapidly than antiprotons. In heavy- 
ion collisions, the p/p ratio decreases with system size, 
and falls more rapidly at mid-rapidity than at backward 
rapidity. This is exactly what is expected from increased 
stopping power, and therefore larger net baryon density, 
with larger colliding system size. 



VII. DISCUSSION 

In heavy-ion collisions, we use the notions of chemical 
freeze-out and thermal freeze-out to describe the moment 
when the inelastic and elastic interactions stop. In this 
section, we will describe the parameters extracted from 
the measured spectra and yields. 
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FIG. 8: Collision system dependence of charged pion, charged kaon, proton, and antiproton dN/dy (top figures), and ratio of 
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antiproton dN/dy are corrected for feed-down. The values of dN/dy in p+A and S+A collisions are from Ref. [|. Statistical 
errors and systematic errors are shown as vertical and horizontal bars above and below each point, respectively. 



A. Transverse momentum distributions 

The transverse momentum distributions of hadrons 
have been measured by several experimental groups at 
BNL-AGS and CERN-SPS including NA44 |j. The 
rriT distributions are studied from the viewpoint of collec- 
tive flow, with local thermal equilibrium at mid-rapidity, 
based on Refs. |]| [l]|. The NA44 collaboration has 
reported a thermal freeze-out temperature T t h ~ 140 
MeV fU in Pb+Pb collisions at v^vw ~ 18 GeV. Re- 
cently two-pion interferometry results have been com- 
bined with single particle spectra to measure the trans- 
verse flow and thermal freeze-out temperature |i~4| 20 . 
These studies find T th « 120 MeV |S] and 100 MeV [§0|. 
Evaluating those reports equally, one concludes that the 
thermal freeze-out temperature is in the range 100-140 
MeV. 



B. Particle yields 

If we assume a thermodynamic description is valid for 
these systems, then at chemical freeze-out, the density of 
different particles can be characterized by macroscopic 
parameters: the chemical freeze-out temperature and 
chemical potentials. Here, we adopt a chemical freeze- 
out model based on Ref. H|, using the grand canoni- 
cal ensemble. We assume no difference between u and 
d quark chemical potentials. Therefore, the hadron gas 



is described by a chemical freeze-out temperature (T c h), 
light (u and d) quark potential ((J>q), strange quark po- 
tential (/i s ), and strangeness saturation factor (7 S ). The 
density of a particle i in the hadron gas is given by: 



Pi = 7 ; 



2tt 2 ' 



T, 



ch 



(2) 

where m$ is the mass of the hadron i, gi is the number 
of spin-isospin degree-of-freedom, K 2 is the second-order 
modified Bessel function and, 

X q = exp(/j, q /T ch ), A s = exp(/i s /T ofc ). 

Note that we use the Boltzmann approximation for all 
hadrons except pions, where the Bose distribution is ap- 
plied. The potential [i q is for u/d/u/d quarks, and /i s is 
for s/s quarks. Qi and Si correspond to the net number 
of valence u/d quarks (Qi = (u — u + d — <i) i ), and s 
quark (sj = (s — s)J of particle species i, respectively. 
The factor j s (0<7 S <1) is introduced to allow of pos- 
sibility of incomplete chemical equilibration for strange 
particles 0|. 7 S =1 signifies full strangeness chemical 
equilibrium. The power factor of 7 S , (s + s),, is total 
number of s and s quarks in the particle. Ideally the sys- 
tem would live long enough that strange quarks would 
become fully equilibrated with the light quarks and 7 S 
would be unity. However, calculations in Refs. |^3|, |24|, |2^] 
support 7 S <1. Therefore, we adopt j s as a free param- 
eter. The particle densities are computed by Eq. (0) for 
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FIG. 9: Collision system dependence of the chemical freeze- 
out parameters. 



TABLE III: The macroscopic parameters for the description 
of the chemical freeze-out near central rapidity for Pb+Pb, 
S+A, and p+A collisions from the NA44 data. The values of 
dN/dy from Refs. [^, El are used to obtain the parameters in 
the p+A and S+A collisions. 





T c h 


fJ,q 


/'* 








[MeV] 


[MeV] 


[MeV] 


7« 


X 2 /DOF 


Pb+Pb 


141+5 


63+7 


19+9 


0.74+0.05 


0.72/1 


S+Pb 


139+6 


58+6 


19+15 


0.61+0.09 


2.6xl0" 4 /l 


s+s 


138+4 


54+4 


20+11 


0.63+0.06 


3.8xl0- 5 /i 


p+Pb 


133+4 


30+4 


10+12 


0.46+0.05 


l.OxlO -8 / 1 


P+S 


137+5 


31+4 


21+11 


0.58+0.07 


6.0xl0~ s /l 


p+Be 


138+5 


26+4 


8+12 


0.49+0.06 


3.9xl0~ 4 /l 



the hadron gas including known resonances up to mass 
of 1.7 GeV/c 2 , and the effect of broad resonance mass 
widths is taken into account pl| . The resonance decay 
to lower mass hadrons after chemical freeze-out is also 
taken into account. Therefore, the final particle densities 
for the pions, kaons, and protons are given by the sum of 
the direct production and the decays of the resonances. 
Here, we apply the above procedure for the particle yields 
observed near mid-rapidity by NA44 for p+A, S+A, and 
Pb+Pb collisions |g || . The model is fit to the data using 
MINUIT Jp). 

Figure gand Table III shows the collision system de- 
pendence of the parameters. The chemical freeze-out 
temperatures (T c h) are approximately 140 MeV in p+Be 
through Pb+Pb collisions: Within errors, there is no de- 



pendence on system size. C omp aring to other results, the 
T ch is consistent with Ref. (25) (7^=133-144 MeV) and 
smaller than Ref. || (7^=158+3 MeV) and Ref. Q 
(T cfc =168±2.4 MeV). As reported in Refs. |lj |§ there 
is a rapidity dependence of the chemical freeze-out pa- 
rameters. Those calculations show that a lower value of 
the temperature is obtained when only mid-rapidity data 



are used. The limited number of particle species used in 
this analysis may also influence the value of T c h . 

The values of n q increase with the size of the collision 
system. Since the chemical freeze-out temperature is in- 
dependent of collision systems (p+A and A+A), the sys- 
tem in p+A collisions has smaller fi q /T c h, that is, smaller 
baryon density than A+A collisions. fi s is 10—20 MeV 
for p+A through Pb+Pb. /z s should be zero in a equili- 
brated quark gluon plasma and on the phase boundary 
of a hadron gas |p2[ p9| |3C[| . Consequently the small 
value we obtain may suggest that chemical freeze-out 
takes place near the phase transition. 

Finally we find that 7 S increases logarithmically with 
dN c h/dy from 0.49 ± 0.06 in p+Be collisions to 0.74 ± 
0.05 in Pb+Pb. Thus, we seem to be approaching full 
chemical equilibrium for strange particles in heavy-ion 
collisions. The use of the grand canonical ensemble for 
analysis of strange particle production in p+A collisions 
does not take into account phase space effects arising 
from production of s and s pairs. Consequently, 7 S is 
expected to be lower than in heavy-ion collisions where 
more strange quark pairs are produced |31, p2| . 



VIII. CONCLUSION 

The NA44 experiment has measured single particle dis- 
tributions for charged pions, kaons, and protons as func- 
tions of transverse mass (mr) near mid-rapidity in 158 
A GeV/c Pb+Pb collisions. We studied particle ratios 
near mid-rapidity from the viewpoint of chemical freeze- 
out for p+A (A = Be, S, and Pb), S+A (A = S, Pb), 
and Pb+Pb collisions. We find that chemical freeze-out 
occurs at a temperature of «140 MeV for all collisions. 
The values of the light quark chemical potential (fj, q ) are 
approximately 30 MeV in p+A and 54 to 63 MeV in 
S+S to Pb+Pb collisions. The trend of \i q shows that 
the baryon stopping power increases with collision system 
size. The values of the strange quark chemical potential 
are about 10 to 20 MeV. The strangeness saturation fac- 
tor is 0.61 to 0.75 in heavy-ion collisions and 0.46 to 0.58 
in p+A, indicating that p+A collisions are further away 
from strangeness equilibration than A+A collisions. 
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